Role of elastic scattering in electron dynamics at ordered alkali overlayers on Cu(lll) 
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Scanning tunneling spectroscopy of p(2 x 2) Cs and Na ordered overlayers on Cu(lll) reveals 
similar line widths of quasi two-dimensional quantum well states despite largely different binding 
energies. Detailed calculations show that 50 % of the line widths are due to electron-phonon scat- 
tering while inelastic electron-electron scattering is negligible. A frequently ignored mechanism for 
ordered structures, i.e., enhanced elastic scattering due to Brillouin zone back folding, contributes 
the remaining width. 
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In many electron systems ubiquitous interactions be- 
tween quasiparticles such as electron-electron, electron- 
phonon or electron-magnon scattering are present and 
cause inelastic scattering. A fingerprint of this dynamics, 
which occurs on a femtosecond timescale, is the width of 
spectroscopic lines. Often, elastic scattering is significant 
too, except for well defined surface states at special points 
in the surface Brillouin zone (SBZ). Thus, surface reso- 
nances exhibit a finite width owing to a degree of mixing 
with the bulk continuum. Defect scattering can also in- 
crease the line width of electron states at artificial atomic 
arrays or at electron confining island Adsor- 
bates on a surface lead to additional scattering channels 
that may change the dynamics of electrons at surfaces 
and at the same time can be an important source of elas- 
tic scattering for electrons |5j. An intriguing situation 
occurs when an extended ordered superstructure changes 
the periodicity of a surface, as in the case of submono- 
layer coverages of alkali atoms on (111) noble metal sur- 
faces. The SBZs of the substrate and of the overlayer 
have different sizes. The resulting back folding of the 
bands can lead to large effects in the dynamics of elec- 
trons at the interface, like the appearance of an elastic 
width of quantum well states (QWS) when energy gaps 
are closed in certain regions of the SBZ. 

Inverse photoemission and two-photon photoemission 
are widely used techniques to probe the dynamics of hot 
electrons at surfaces |6j, |j| ■ Recently, the different mecha- 
nisms that contribute to the elastic and inelastic scatter- 
ing of hot electrons at Cu(001) surfaces with a low cov- 
erage of Cu adatoms have been identified using these 
techniques. The conclusions are deduced from informa- 
tion taken over the whole surface area without probing 
the local electronic structure. Scanning tunneling spec- 
troscopy (STS) resolves the local electronic structure of 
surface areas free from defects. It also can provide de- 
tailed information about both occupied and unoccupied 



states. The lifetime broadening of surface states can be 
extracted from STS measurements by doing a line shape 
analysis Wi JToL fllL or from electron standing wave 
patterns |13l Il4. Il5| . The progress of these experimental 
approaches clearly requires refined theoretical methods if 
one aims to explain experimental data at a quantitative 
level. 

Here, in a combined experimental and theoretical 
study, we analyze the processes that contribute to the 
width of QWSs at Cs and Na overlayers on Cu(lll). We 
achieve quantitative agreement between low temperature 
STS experiments and simulations of tunneling spectra us- 
ing first principles calculations of the electronic structure 
of the surface and of the line width. The width of the 
QWS is shown to have significant contributions from the 
elastic width that appears as a result of the coupling to 
substrate states and from the inelastic electron-phonon 
scattering whereas inelastic electron-electron interaction 
is of minor importance, contrary to the situation ob- 
served in studies of the hole dynamics at the (111) surface 
of noble metals 01 • 

Clean Cu(lll) surfaces were prepared by standard Ar + 
bombardment and annealing cycles. Cs and Na were de- 
posited from commercial dispensers [3] onto the sam- 
ple at room temperature at rates of « O.SMLmin -1 as 
monitored with a quartz microbalance. We define one 
monolayer (ML) as one alkali atom per one Cu atom of 
the unreconstructed Cu(lll) substrate. In the Cs case, 
a clear p(2 x 2) low-energy electron diffraction (LEED) 
pattern at room temperature [ItI was used as ad- 
ditional indication of a 0.25 ML coverage. After prepa- 
ration, samples were transfered to a scanning tunneling 
microscope (STM) and cooled to temperatures of 4.6 — 9 
K. Spectroscopy of the differential conductance (dl/dV), 
was performed on homogeneous areas using a lock-in am- 
plifier and adding a sinusoidal modulation (6 — 14mV pp , 
10 kHz) to the sample voltage. 
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Figure ^ displays d//dV-spectra (dots) of p(2 x 2) 
Cs and Na layers. Sharp rises of the conductance at 
Eq = 40 meV and E = 410 meV correspond to the 
minima of the quantum well state bands of Cs and Na 
overlayers, respectively. The width of the rise, which is 
related to the inverse lifetimes of the these states is 
similar in both cases. This observation is surprising in 
view of the large difference in binding energy which, in 
a simple picture, should cause widely different electron- 
electron (e — e) scattering. However, this appears not to 
be the case. Our calculations, whose results agree well 
with the experimental data (solid lines in Fig. 0), show 
that inelastic e — e scattering is unimportant and that 
electron-phonon (e — ph) scattering explains only 50 % 
of the line widths. Enhanced elastic scattering due to 
Brillouin zone back folding, turns out to contribute the 
missing 50 % of the line width. 

The surface electronic structure of both systems 
has been calculated using the density functional code 
VASP 19]. Projected Augmented Waves [2(| are used 
to describe the ionic cores and a plane wave basis set 
for the valence electrons. The exchange and correlation 
potential is treated in the generalized gradient approx- 
imation [2]]]. We have checked the convergency of the 
calculations with respect to the thickness of the slab up 
to 10 layers of Cu. The fc-point sampling of the irre- 
ducible wedge of the surface Brillouin zone of the (2x2) 
superstructure is done using 14 special points for the 
self-consistent calculation and later on we use up to 154 
points to achieve high enough resolution of the bands. 
These slab calculations provide the structure for the re- 
cursive calculation of the elastic width in the semi-infinite 
medium. The tunneling conductance is calculated with 
a model that assumes a smooth density of states of the 
tip around the Fermi level. Then, at small positive sam- 
ple bias, the energy dependence of the conductance is 
determined by the sample density of states 1221 . The de- 
tails of the model are described elsewhere [23J . In brief, 
we introduce an energy width (T m /2) of the states when 
calculating the imaginary part of the surface Green func- 
tion, i.e., the surface local density of states (LDOS). This 
width corresponds to a Lorentzian weight for each state 
such that the surface LDOS is given by: 

p( e .^P*^ Vj"+ 2 (W2)' (1 » 

The eigenvalues (e m ) and eigenfunctions (\& m ) are 
taken from VASP, while the energy width (r m /2) is de- 
termined as described below. The sum over m includes 
both a summation over the band index and an integra- 
tion over parallel momentum. The LDOS is calculated 
at 3 A from the alkali overlayer, where the electronic cor- 
rugation is negligible, and it is dominated by the QWS 
band in the relevant energy range. A /c-point interpo- 
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FIG. 1: Comparison between measured d//dV vs. V curves 
(dots) for p(2 x 2) structures of (a) Cs and (b) Na on Cu(lll) 
along with calculated results (lines). The solid lines include 
both the elastic and inelastic contributions to the energy 
width, while the dashed lines include only the inelastic one. 



lation of the QWS is used to achieve the required meV 
resolution of the band [i^j ]. 

For comparing the measured d//dV data for the p(2 x 
2) Cs structure on Cu(lll) and the simulated spectrum 
we return to Fig.^,. No fitting parameter has been used 
apart from a scaling factor for the calculated conduc- 
tance. The agreement between the data points and the 
solid curve clearly shows the importance of the elastic 
width in determining the shape of the spectrum at the 
onset around 40 meV. Fig. \\]p displays a similar compar- 
ison for a p(2 x 2) Na film on Cu(lll). In this case, the 
calculated curves have been shifted by +18 meV to match 
the experimental binding energy. This small shift is con- 
sistent with the expected accuracy of density functional 
calculations [24J. The shape of the curve at the onset 
also agrees well with the measured data in this case. 

The adsorption of alkali adlayers in p(2 x 2) super- 
structures leads to a modification of the projected bulk 
band structure of the substrate as shown in Fig. [21 along 
the r — M direction. The inset displays the SBZs for a 
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p(2 x 2) overlayer and for a (1 x 1) substrate. In the 
over layer case the energy gap near Ep of clean Cu(lll) 
surface is closed due to the folding of bands originally 
situated along the M — M Cu and M Cu — K' symmetry 
directions. The corresponding extra corrugation of the 
Cu(lll) surface due to the coupling between the p(2 x 2) 
adsorbate layer and underlying Cu bulk states causes a 
degree of mixing between the alkali induced QWS at T 
and substrate Cu states at the Mc u point of the (lxl) 
Brillouin zone. This serves as a source of elastic scatter- 
ing, which requires a quantitative analysis. 

To calculate the elastic width of QWSs we compute 
the Green function of the combined system (adlayer plus 
semi-infinite substrate) with high precision in the surface 
region. Here, we use the SIESTA code l25l w hich utilizes 
a basis set of numerical atomic orbitals |26( . The results 
for the electronic structure are almost identical to those 
obtained with VASP. With this local basis set the ele- 
ments of the Hamiltonian and overlap matrices between 
atoms that are far apart (beyond «15 A in our case) 
are strictly zero. The infinite system can now be divided 
in groups of layers (" slices" ) that only interact with the 
nearest neighbor groups. The Hamiltonian matrix ele- 
ments (and overlaps) for orbitals in the surface slice, and 
its interaction with orbitals in underlying layers, are ob- 
tained from the slab calculation. The interactions in the 
inner slices are taken from a calculation of bulk Cu. A 
common energy reference is set by aligning the Fermi 
levels of both calculations. We now use the recursive 
relation 

G«(o; ) k||)[^(k||)-^ fc (k||)]=4 (2) 

to obtain the Green function in the surface region. We 
use a method similar to that described in Ref. j2?| to 
iteratively solve this equation. The matrix elements are 
defined such that G(r, r'; w, kii)= (a;, k|| ) fa(r)(j)j(r'), 
where the fa are the atomic orbitals and the sum over re- 
peated indices is assumed. Our strategy is then to project 
the Green function onto a wave packet \& r localized in 
the surface region. The function Im{< ^r\G(lo)\^r, >} 
exhibits well defined peaks at energies cr corresponding 
to the QWS of the overlayer. The shape and position of 
these peaks are independent of the election of 4" r, as far 
as ^r has predominant contributions from the orbitals 
in the adsorbates. The peaks have a Lorentzian shape 
with a width and position that can be easily extracted 
using a least square fit |28j ] . Carrying out this procedure 
for the two systems studied here we find values of the 
elastic width at T equal to 9.4 meV and 7.4 meV for Cs 
and Na systems, respectively. 

For both systems of interest we estimate the inelastic 
electron-electron (e — e) contribution (r e _ e ) to the QWS 
linewidth by using a model in which the one-electron po- 
tential is constant in a plane parallel to the surface and 
varies only in the direction z perpendicular to the sur- 
face. We use here the potential proposed in Ref. [2j| 
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FIG. 2: Dispersion of p(2 x 2) Na and Cs quantum well states 
in the T — M direction. Dark grey indicates bulk Cu states of 
the (1 x 1) SBZ, while light grey shows the closing of the gap 
around F after back folding from the M — Mcu and Mcu — K' 
symmetry directions defined in the inset. 

for thin films adsorbed on a Cu(lll) substrate with pa- 
rameters fitted to reproduce the binding energy of the 
QWS in the Cs (42 meV) and Na (408 meV) adlayer on 
Cu(lll). r e _ e is calculated within the GW approxima- 
tion by employing eigenfunctions and eigenenergies 
obtained with this model potential. The computed T e _ e 
values of less than 0.1 meV and 0.4 meV for a QWS in 
the Cs and Na systems, respectively, demonstrate that 
inelastic e — e scattering provides only a small contribu- 
tion to the QWS linewidth. This latter fact reflects the 
difference between electrons and holes: for electrons at 
the bottom of the QWS the intraband channel is com- 
pletely absent, while it is most important for holes [Icj . 

Another inelastic contribution, the electron-phonon 
(e — ph) one, is found to be important for the stud- 
ied systems. While there is no detailed information on 
phonon dispersion of these systems, one can refer to he- 
lium atom scattering measurements for a monolayer of 
Cs on Cu(100) 30]. Assuming that this phonon struc- 
ture is also valid for Cs on Cu(lll) we evaluate V e -ph by 
using the Einstein model [3JJ to treat a flat phonon mode 
with an energ y o f 7.5 meV around the Y point and the 
Debye model |3l| for the two acoustic modes (Rayleigh 
and longitudinal ones) with maximum energy of 6 meV. 
The e — ph coupling parameter A = 0.18 is taken from 
Ref. j3]| as the maximum bulk value for A. With this 
we estimate T e _ p f l = 7.5 meV. For Na on Cu we also 
assume that the measured maximum phonon frequency 
for 0.44 ML of Na on Cu(100) can be applied to p(2 x 2) 
Na/Cu(lll). Following Ref. ^ we use A = 0.24 and 
Debye energy Tuvd = 18 meV yielding r e _ p /i = 9.0 meV. 

An overview of the different contributions to the line 
width of the QWS at the T point is presented in Table 
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^elastic 


Te-e 




-T total 


Cs 


9.4 


<0.1 


7.5 


17 


Na 


7.4 


0.4 


9 


16.8 



TABLE I: Line width contributions (in meV) of various scat- 
tering channels. F e i astic elastic, F e _ e electron-electron, and 
r e _ p h electron-phonon scattering. Y to tai is the resulting total 
width. 



The elastic contribution to the energy width accounts 
for approximately 50 % of the total width. The inelastic 
e — ph contribution is comparable. However, the inelastic 
e — e contribution is negligible. 

In summary, our combined experimental and theoreti- 
cal study shows that the dynamics of electrons at ordered 
overlayers is determined by the elastic width that appears 
due to the closing of gaps after back folding the bands 
and by inelastic e-ph scattering. However, inelastic e-e 
scattering does not play a role. The relative importance 
of the elastic contribution for other ordered superstruc- 
tures will depend on details like the density of states at 
given points of the substrate SBZ. Any change of period- 
icity on going from the surface to the overlayer structure 
will lead to a modification of the scattering channels that 
determine the dynamics of electrons at the interface. 
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